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Abstract 

We study the nuclear structure function i<2^ and its logarithmic derivative in the 
high energy hmit (small x region) using the Color Glass Condensate formalism. In 
this limit the structure function F2 depends on the quark anti-quark dipole-target 
scattering cross section NF{xbj,rt, ht). The same dipole cross section appears in single 
hadron and hadron-photon production cross sections in the forward rapidity region 
in deuteron (proton)-nucleus collisions at high energy, i.e. at RHIC and LHC. We use 
a parameterization of the dipole cross section, which has successfully been used to 
describe the deuteron-gold data at RHIC, to compute the nuclear structure function 

and its log derivative (which is related to gluon distribution function in the 
double log limit). We provide a quantitative estimate of the nuclear shadowing of 

and the gluon distribution function in the kinematic region relevant to a future 
Electron-Ion Collider. 



1 Introduction 



The Color Glass Condensate (CGC) formalism [Ij has been quite successful in describing 
aspects of particle production in high energy hadronic/nuclear collisions. It has been 
used to describe data from lepton-proton (nucleus) Deeply Inelastic Scattering (DIS) to 
hadron multiplicities and transverse momentum spectra in nucleus-nucleus and deuteron- 
nucleus collisions [2l|. The recently observed suppression of single hadron pt spectra in the 
forward rapidity region in deuteron-nucleus collisions at RHIC [3j was predicted by the 
formalism [1]. Since then, quantitative results which describe the pt spectrum of single 
hadron production in dA collisions have become available [5], [6], [7] . Therefore it is desirable 
to apply the knowledge gained from dA collisions at RHIC to other processes in order 
to further clarify/constrain the region of applicability of CGC formalism to high energy 
processes [8|. 

In this work we consider nuclear shadowing of structure function and its log 
derivative using the results of our previous work on pion production in dA collisions. In 
[HI E] a hybrid approach to dA collisions was used to investigate pion production in dA 
collisions (for an alternative approach which treats both the projectile and target by CGC 
formalism, see |10]). The main result of that analysis |6] can be summarized roughly as 
the following: in the mid rapidity region, as one goes from low to high pt, one goes from 
the saturation region through the scaling region to the pQCD (DGLAP) region. The 
applicability of CGC in mid rapidity region of RHIC is limited to less than a few GeV in 
transverse momentum above which one is in the large x region. It is worth mentioning 
that the qualitative conclusion reached in [7] is different from those of |i6j due to a different 
parameterization of the dipole cross section which enters the particle production cross 
section. However, the dipole model in [7] seems to have an unphysical dependence on 
transverse momentum at high pt. Therefore, we will use the parameterization of the dipole 
cross section in p] (referred to as the DHJ model) in this work. 

In the most forward rapidity region at RHIC (and LHC), one is in the (target) small 
X kinematics {pt is limited by kinematics) but in the projectile large x region where the 
quarks in the projectile contribute more than the projectile gluons (except possibly at very 
low Pt where they are comparable). Therefore any description of the forward rapidity data 
must include the effects of the projectile quarks scattering on the gluon field of the target 
nucleus (or proton). 

Using the DHJ model of the dipole cross section employed in [6j, we calculate nuclear 
structure functions and their shadowing. The main ingredient in the structure function F2 
is the (quark anti-quark) dipole cross section which also appears in single hadron production 
cross section in dA collisions. Specifically, we investigate how the structure function 
computed using the DHJ model of the dipole cross section describes the nuclear shadowing 
data as measured by the NMC collaboration [11]. We then make predictions for in 
a broader kinematic region which may be covered by the proposed Electron-Ion Collider 
(EIC). We then consider shadowing of the nuclear gluon distribution function. In a high 
gluon density environment the gluon distribution function (defined as the leading twist 
matrix element of gluon field operators) is not what goes into a physical cross section. 
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rather, it is n-point functions of Wilson lines. Therefore, we consider the log derivative 
of the structure function which in the high (double log) limit is proportional to 
the gluon distribution function. We then compare the ratio of these log derivatives, for a 
nucleus and a proton, with the results available in the literature. 

2 The structure function 

In the small x limit, the structure function can be written as a convolution of the 
probability for splitting of a virtual photon into a quark anti-quark pair with the probability 
for the scattering of the quark anti-quark pair on the target so that [12] 

F2 = aQ^ J drtn J dzNp{x,rt) {^fi{z) K^,{ert) + fo{z) K^ien)^ (1) 

where a = is an overall constant and fi{z) = z{l — z)[z'^ + (1 — z)"^], fo{z) = 

4 — z)^. The sum is over quark flavors (we are taking three flavors of massless quarks) 
and we have defined = z{l — z)Q'^. All the QCD dynamics in the small x limit is encoded 
in the dipole cross section Np{x,rt) which is also the building block of the single hadron 
production cross section in the forward rapidity region and satisfies the JIMWLK/BK 
evolution equations. 

It would be ideal to solve the JIMWLK/BK equations in order to obtain the form of 
the dipole cross section Np and use it in eq. ([T]) (see [13] for a recent analysis of the proton 
structure function using a solution of the BK equation with running coupling). However, 
the numerical solutions of JIMWLK/BK equations are known to be quite sensitive to the 
initial condition in the evolution equation unless one is at extremely small x. It is also 
expected to receive large corrections from higher order corrections. Therefore it is more 
practical to parameterize the dipole cross section in a form which captures the essential 
dynamics of the saturation physics. This approach has been employed before to investigate 
structure functions at HERA [14j and particle production at RHIC [5| [Tf [T5]. 

In this work we use the parameterization proposed in [6J which was fit to the single 
hadron production data in dA collisions at RHIC at rapidity y = 3.2. It was then used 
to predict the data at rapidities y = and y = 4. For a successful description of the 
mid rapidity data one had to exclude the contribution from x > 0.05. This is consistent 
with the lessons learned from HERA where CGC motivated models can not fit the data 
(for proton targets) for x > 0.01. Therefore, we will restrict ourselves to the region where 
X < 0.05 which puts a severe limit on the number of data points available on nuclear 
shadowing at reasonably large (we consider > 1 GeV^ only). There was also a need 
for a K factor in description of the single hadron spectra in dA collisions which was large 
(in mid rapidity) but pt independent (see [16] for observables which are independent of 
this K factor). This is usually attributed to higher order corrections which are reduced in 
magnitude as one goes to the forward rapidity region as was the case in [6]. Here we do 
not expect a K factor since we are dealing with a fully inclusive quantity. Furthermore, we 
do not include the data points measured the E665 collaboration since there seems to be 
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a discrepancy between the NMC and E665 shadowing results when comparing structure 
functions for large nuclei (such as gold or lead) to that of deuteron which may be due to 
the systematic errors of the E665 measurement [T7] . 



2.1 The dipole cross section and shadowing of FA 
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In DHJ parameterization [6j, the fundamental (quark anti-quark) dipole scattering cross 
section is given by (an overall factor of (Tq = 23 mb is factored out in this expression but 
included in the numerical results) 

Nj.{x,n) = 1 - exp [-i[r2Q2(^)]7(-,n)j (2) 

where the anomalous dimension 7 is given by (for details see [6]) 

7irt,y) = 7s + A7(r,?/) 

where y = logl/x and d = 1.2 (see [6J for details of the parameterization). Due to 
the anomalous dimension 7, there is a strong "leading twist" shadowing of the target 
gluons encoded in this formalism [18]. Using this, we compute the minimum bias structure 
functions for a proton and a nucleus (taken to be gold) in the kinematic region covered by 
the NMC experiment while keeping x < 0.05 and > 1 GeV"^. 

Our result for nuclear shadowing is shown in Fig. ([T]). The lower (solid) line is the 
ratio of F2 / F'2 where both structure functions for proton and gold are calculated using 
this formalism. The upper (dashed) line is the ratio F^jF^ where we have included an 
overall 10% uncertainty to account for the effect of deuteron shadowing (about 2 — 3% in 
this kinematics [H]), the fact that the experimental data is for lead rather than gold and 
the fact that our calculation of the proton structure function overestimates the HERA data 
by about 10% [U]. All these effects would make the calculated shadowing ratio increase. 
The fact that we are over estimating the proton structure function is not unreasonable 
since the DHJ parameterization was developed and tested for nuclear rather than proton 
targets. We emphasize that there is no new parameters introduced for the computation of 
the structure function. Alternatively, we could have divided the nuclear structure function 
F2 calculated using this formalism by the proton structure function measured at HERA. 
Since we will consider the shadowing of the "gluon distribution function" defined as the 
log derivative of F2, we have decided to divide by the calculated value of rather 
than the measured value and thus have included the upper line as a more realistic estimate 
of the shadowing effect. As is clear, the difference between the data points is largest (of 
order of 15%) at the highest x considered [x = 0.05) and gets much smaller (a few percent) 
toward the lowest x in the figure {x = 0.01). 

In Figure ([2]) we show the nuclear structure function F.f as a function of at different 
values of x. We show our results for values of up to 400 GeV^ which was an estimate for 
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the absolute upper limit of the scaling region at HERA [20] . The analysis in [B] shows that 
the scaling region may actually be smaller. Furthermore, the calculation of is more 
reliable in our case since we don't have to consider the proton structure function which our 
calculation overestimates bv about 10%. 




0.05 



Figure 1: Shadowing of structure function compared to NMC measurements. 



In Figures ([3]) and (jlj) we show the shadowing effect for the structure function at 
different kinematics. In both figures we have plotted the ratio of F2 for minimum bias 
gold nucleus and that of a proton and have ignored the shadowing of deuteron structure 
function at small x (which may eventually be about 10%) and the fact that our calculation 
overestimates the proton structure function at higher by about 10%. Therefore, the 
lines in both figures should be understood to be a lower limit which could be pushed up 
by 10 — 15% in a more detailed study. 

It is worth noting that one could fit the dependence of the lines shown in Fig. ([3]) 
by a function of the form cq + ciLog / Ql with constants cq, ci being the same (with 2% 
accuracy) for all three lines shown (the deviation is largest at smallest which may be 
a sign of reaching the saturation region). Whether the saturation scale appearing in the 
log is that of a proton or nucleus can not be determined from the fit since the constant Ci 
is much smaller than cq. It is tempting to take this as a sign that the shadowing in this 
kinematics is "leading twist", consistent with the estimates of the scaling region where the 
dynamics is that of BFKL with saturation boundary (we note that the saturation scale of 
a minimum bias gold nucleus is 1.6 GeV at x = 3 x 10~^ and IGeV for a proton). On 
the other hand, in the much larger kinematics shown in Fig. (jlj) this functional form does 
not provide a very accurate fit which may be a sign that we are outside the scaling region. 
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Figure 2: Nuclear structure function Fi- 
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Figure 3: Nuclear shadowing of the structure function F^. 
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Figure 4: Same as in Fig. ([3]) but in a larger kinematic region. 



It would be very interesting to make analytic estimates for the extent of the scaling 
region in this case, however, this is a bit complicated due to the z integration in eq. ([T]). 
In principle, the value of photon virtuality determines whether one is in the saturation 
region Q -C Q^, in the scaling region Qs <^ Q Qll ^qcd or in the pQCD region where 
Q1/Aqcd ^ Q- The the integration over the dipole size can be broken up into 4 distinct 
regions 



The first region is that of pQCD where rtQ <^ 1 and VtQg ^ 1. Then Np ~ f^tQ's 
K\ ~ and Kq is subleading. In this case the rj and z integrations can be performed 

analytically and are finite. A similar thing happens in the saturation region (the last piece 
in eq. (jlj)). In this case, rtQ ^ 1 and rtQs^ ^ and one can approximate Np ~ 1 and 
Kq ~ Kl ~ {1/ ert) exp{—2ertQ). However, in the other two regions, the z integration 
diverges as 2; — > 0, 1 and is not under control (having massive quarks would regulate this 
divergence). One can still perform the integrals analytically (for fixed 7) but its general 
form is not very illuminating. 

2.2 Shadowing of gluons 

In this section we consider nuclear shadowing of the gluons. In pQCD and at the leading 
twist level, the gluon distribution function xG is proportional to the log derivative of 
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F2 structure function, 



xG{x,Q') ^ jj^^F,{x,Q') (5) 

This relation can be used to extract the gluon distribution function experimentally in DIS. 
The Color Glass Condensate formalism, however, is not based on a twist expansion and 
the relevant degrees of freedom (which appear in physical cross sections) are not gluons 
but Wilson lines, path ordered exponentials of gluon fields which take multiple scattering 
on the target into account. In a typical cross section computed in the CGC formalism, one 
encounters n-point functions of Wilson lines (fundamental or adjoint), where in case of the 
structure function F2, only the two point function of (fundamental) Wilson lines appears. 
Therefore, in order to provide a comparison with results obtained from other approaches 
based on leading twist coUinear factorization (where shadowing is put in by hand in the 
initial condition), we consider the log derivative of the F2 structure function and take 
that to be a measure of shadowing of gluons in the standard language. 

Differentiating the structure function F2 as given in eq. ([1]) with respect to log Q^, 
we get two terms; one is coming from differentiating the overall Q'^ term and the second 
contribution is coming from the dependence of the Bessel functions. We get 



dFo 



aQ' J dnn J dzNF{x,n)[2fo{z)K^{en) - 



d log Q"^ 

ert [fo{z) + fi{z)] Ko{en)Men) + fi{z) K^en) }■ (6) 



Using the expression for Np given by ([2]) we can now calculate the log derivative of the 
structure function. In Figs. ([5l [6]) we show R^g defined as 

dF^ dF^ 
" dlogQ^' dlogQ^' ^ ' 

Here we have ignored shadowing of gluons in a deuteron. We are also using our formula 
to calculate this log derivative for a proton target rather than using xg from available 
parameterizations such as CTEQ or MRST for self consistency. Because of this we may 
underestimate the amount of shadowing by ~ 10 — 15% in the small x but large region. 
This would mean, for example in Fig. ([5]) that this ratio is very close to 1 at the highest 

shown. At the largest x considered [x = 0.01), the amount of gluon shadowing is 
about 10 - 15% at ~ lOGeV^ and 25% at the = IGeV^. It goes away by of 
several hundred GeV^. One should keep in mind that the shadowing results shown are 
most likely an overestimate (the real ratios are a bit larger) since first, we are ignoring 
deuteron shadowing and second, we are overestimating the gluon distribution in a proton 
by possibly 10%. At smaller x and Q"^ = 1 GeV^ gluon shadowing is about 30% and again 
goes away at higher Q'^. The Q"^ dependence of shadowing at fixed x is shown in Fig. ([H]). 
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Figure 5: Shadowing of the gluon distribution function as defined in eq. ([7]). 




Figure 6: Same as in Fig. ([5]) but vs. x. 
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3 Summary 



We have used the parameterization of the quark anti-quark dipole successfully employed 
in [6] to investigate the predictions of the model for shadowing of minimum bias nuclear 
structure functions. We have compared our results to the available data from NMC col- 
laboration in the smallest x region accessed in that fixed target experiment while keeping 
the values of larger than IGeV"^. Our results are within a few percent of the data 
at the smallest x bin (~ 0.0125) and are off by about 15% at the highest x considered 
(~ 0.045). Defining gluon distribution function as the log derivative of the F2 structure 
function, we calculated the shadowing of gluons and have made predictions for its x and 
Q"^ dependence which could be tested in a future lepton- nucleus collider experiment. 

It is interesting to compare our results with those from other approaches. Typically 
there are two distinct parton-based approaches (see also [21] for an alternative approach) 
to particle production in high energy collisions and nuclear shadowing of structure func- 
tions. First is based on the leading twist coUinear factorization and DGLAP evolution of 
distribution functions. In this approaches one has to put in shadowing by hand as a modi- 
fication of the distribution functions at some scale Qq which are then evolved according to 
DGLAP equations [221 1231 [23|. Attempts to simultaneously fit the RHIC data on forward 
rapidity dA collisions and the nuclear structure functions seem to lead to an almost 100% 
shadowing of gluons in a gold or lead nucleus at small x which seems quite unnatural [23] 
(this ratio seems to be higher in the most recent version EPS09 [21] which uses a NLO 
analysis). Comparing our results for shadowing of gluon distribution function as defined 
in eq. ([6]) to those in [231 [23] , our shadowing ratio is much closer to the rsults of [23] than 
[24j . For example, we predict Rxg ~ 0.65 (for minimum bias gold nucleus) at x = 10~^ and 

~ l.lGeV^ while this ratio in [23] (LO) is ~ 0.5, but ~ in [23] and ~ 0.5 in the most 
recent version in [23] (NLO). The difference in the amount of gluon shadowing between 
[23] and [23] (LO) is quite interesting. The drastic shadowing of gluons in [23] (there are 
almost no gluons in the nucleus) is caused by insisting on describing the RHIC forward 
rapidity data using collinear factorization approach where nuclear shadowing is put in the 
initial gluon distribution by hand. Since having almost no gluons in the nucleus seems 
a bit unnatural, this may indicate that leading twist collinear factorization approach to 
forward rapidity data at RHIC is not valid. The amount of gluon shadowing in the NLO 
version is much less than in the LO version but it is not clear whether the NLO version 
also reproduced the forward rapidity RHIC data. There has also been attempts to include 
the effect of the higher twist corrections to the DGLAP evolution equations [25j but these 
effects would break collinear factorization and make the approach practically useless for 
hadronic/nuclear collisions. 

The second approach is based on the concept of gluon saturation and here also one can 
distinguish between two different approximations to the dipole cross section (see also [26] 
which studied shadowing using scaling ideas). One is motivated by the original approach 
of Golec-Biernat and Wusthoff which is a Glauber like multiple scattering approach to 
shadowing but using the parton language. The shadowing generated here is a "higher 
twist" (suppressed by powers of Q"^) effect. One such work is [13] (last paper) which 
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also investigates the impact parameter dependence of shadowing, but does not include the 
effects of BFKL anomalous dimension 7 and in this sense does not have "leading twist" 
shadowing. A second approach within the saturation picture includes this anomalous 
dimension but using a different model for the dipole cross section ^J. It is also used to 
fit the RHIC data in the forward rapidity data but to the best of our knowledge it has not 
been compared against the NMC data on shadowing of F2 (it has been compared with the 
E665 data however) and shadowing of gluons is not considered. 

The fact that our parameter free calculation of the F2 structure function is so close to 
the experimental values is another indication that the observed suppression of the single 
hadron spectra in the forward rapidity dA collisions is due to gluon saturation dynamics. 
Nevertheless there are some interesting questions which have not been addressed here and 
need further investigation. In this context, the first and most important is probably a 
study of the impact parameter dependence of structure functions including BFKL dynam- 
ics. The DHJ parameterization of the dipole cross section was proposed for minimum bias 
dA collisions. It would be interesting to do a similar analysis for the centrality dependent 
hadron spectra in dA collisions and apply the knowledge to predict the centrality depen- 
dence of the structure functions which could be measured in a future facility. This would 
help clarify the dynamics of shadowing, whether is it power or logarithmically suppressed. 
Since the magnitude of saturation scale is much larger for most central collisions, then the 
saturation region for central collisions extends to higher values of as compared with 
minimum bias collisions where the value of is not that large. This work is in progress 
and will be reported elsewhere [27]. 

One could also differentiate between different approaches (CGC vs. collinear factor- 
ization inspired) to shadowing by considering two particle correlations |i28j. In the CGC 
approach where one describes both the projectile and target using CGC formalism, one 
expects a de-correlation of the two particles as the rapidity separation between the two 
observed particles gets larger. This is due to the small x evolution of the gluon ladder 
between the produced hadrons. In the standard collinear factorized form of particle pro- 
duction where shadowing is put in by hand, there should be no de-correlation. Therefore, 
a detailed study of pt, impact parameter and rapidity dependence of two particle produc- 
tion would go a long way toward clarifying the dynamics of shadowing of nuclear structure 
functions. 
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